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ABSTRACT 

A  satellite-communications  terminal  geolocation  is  analyzed  by  satellite  inter¬ 
ception,  whereby  a  number  of  spaceborne  interceptors  transpond  the  frequency  band 
of  interest  to  a  terrestrial  location  for  processing.  The  general  formulas  for  the  inter¬ 
ception  accuracy  are  summarized,  and  a  strawman  interceptor  system  is  proposed. 
Interception  regions  for  prototypical  terminals  and  satellites  are  calculated  and  the 
results  are  presented  parametrically  as  a  function  of  uplink  SNR.  The  optimum 
angular  separation  of  the  interceptor  satellites  is  found,  and  the  effect  of  nonopti- 
mal  separation  is  discussed.  The  practical  limitations  involved  in  implementing  this 
geolocation  system  are  also  discussed. 
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1.  INTRODUCTION 


This  report  analyzes  an  emitter-location  technique  that  involves  spaceborne  interception  of 
ground-to-satellite  communication  links.  In  this  scheme  a  number  of  interceptor  satellites  transpond 
the  frequency  band  of  interest  to  a  terrestrial  location  for  processing,  as  shown  in  Figure  1.  The 
basic  technique  can  be  used  on  a  variety  of  signal  types:  fixed- frequency,  frequency-hopped,  pseudo¬ 
noise,  etc.  However,  the  focus  throughout  this  report  is  the  theoretical  feasibility  of  geolocating 
terminals  that  are  hopping  in  frequency.  The  technique  is  applicable  at  any  frequency,  but,  to 
simplify  the  formulation,  we  assume  that  the  ground  terminals  have  narrow  antenna  beamwidths 
and  they  operate  in  the  higher  satellite  communication  frequency  bands. 


158125-' 


INTENDED  S/C 


Figure  1.  Components  of  geolocation  system. 


The  time  of  arrival  of  a  signal  having  a  known  structure  can  be  estimated  by  correlating  the 
incoming  signal  with  a  replica  and  finding  the  peak  of  the  correlation  function.  In  the  case  of  a 
signal  having  an  unknown  structure  that  arrives  at  two  separate  receivers,  the  differential  time  of 
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arrival  (DTOA)  can  be  similarly  found  by  forming  the  correlation  function1  of  the  two  receiver 
outputs.  Uncertainty  in  frequency  as  well  as  in  time  leads  to  the  generalization  of  the  correlation 
function  known  as  the  complex  ambiguity  function  [1] 

A(r,f)=  fT  Sl(t)S^(t  +  r)e-^tdt  ,  (1) 

Jo 

where  si  and  S2  are  the  complex  envelopes  of  two  signals  that  have  some  common  features.  |  A  (r,  /)[ 
achieves  its  peak  for  values  of  r  and  /  that  are  the  complements  of  the  relative  time  and  frequency 
offsets  of  the  two  signals. 

Given  the  complex-ambiguity-function  data,  the  use  of  a  reasonable  peak-location  algorithm 
will  yield  estimates  of  the  signals’  DTOA  and  differential  frequency  of  arrival  (also  known  as 
differential  Doppler,  or  DD  for  short).  For  a  reasonably  high  signal- to-noise  ratio  (SNR),  the 
algorithm  will  come  close  to  achieving  the  Cramer-Rao  bounds  for  DTOA  and  DD  estimation  [lj. 
Multiple  hyperboloids  that  correspond  to  constant  values  of  DTOA  or  DD  intersect  the  earth,  each 
forming  a  line-of-position  (LOP).  The  intersection  or  near-intersection  of  multiple  LOPs  yields  the 
geolocation  estimate,  the  accuracy  of  which  can  be  derived  in  terms  of  the  accuracies  of  the  DTOA 
and  the  DD  estimates  [2]. 

A  summary  of  known  formulas  for  the  accuracy  of  the  DTOA-based  estimates  and  the  re¬ 
sulting  geolocation  are  provided  in  the  following  section.  Emphasis  is  placed  on  the  simplifications 
that  ensue  when  the  formulas  are  applied  to  ideal  satellite  interception  of  a  relatively  narrow-beam 
transmitted  signal.2 


frequently,  this  is  also  called  the  cross-correlation  function. 

2Although  similar  results  for  DD-based  geolocation  can  be  readily  developed,  they  will  not  be 
explicitly  considered  in  this  report. 
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2.  CRAMER  RAO  BOUNDS  ON  THE  GEOLOCATION  ERROR 


Let  Bs  be  the  signal  bandwidth,  B  the  noise  bandwidth,  and  T  the  duration  of  the  waveform 
that  is  processed.  Furthermore,  7 ^  is  defined  to  be  the  effective  SNR,  i.e.,  Ps/(NqB),  at  the  output 
of  the  (complex)  cross-correlator.  The  cross-correlator’s  output  SNR  can  be  expressed  in  terms  of 
the  input  SNRs  as 


lout  =  BT 


7i72 

1+71+72 


(2) 


where  71  and  72  are  the  SNRs  at  the  inputs  of  the  two  receivers.  Assuming  the  input  SNRs  are 

A 

equal  (71  =  72  =  7m)-  Equation  (2)  takes  on  the  standard  form  of  the  cross-correlator  input /output 
SNR  characteristic 


lout  =  BT 


o 


1  +  27,n 


(3) 


Making  the  substitution  ^  =  RT7,  Equation  (3)  can  be  rewritten  as 


bt+2(&).„ 


(4) 


This  equation  indicates  the  input /output  SNR  relationship  is  linear  at  high  SNR  and  square-law 
at  low  SNR,  with  a  crossover  at  about  -2  dB. 

Hereafter,  it  will  be  assumed  that  ^  refers  to  the  identical  SNRs  at  the  inputs  of  the  receivers, 
and  that  B  =  1/T  =  Ba  =  the  hopping  frequency.  The  latter  assumption,  which  is  valid  when 
there  is  but  one  chip  per  uplink  hop,3  implies  that  7^  =  (E/No)out\BT=1. 

The  1-ct  error  for  the  DTOA  estimate  is  [1) 


&DTOA 


Vz  1 


3An  example  of  a  one-chip-per-hop  situation  is  frequency-hopped  M-ary  FSK  modulation,  where 
each  hop  consists  of  a  single  FSK  tone. 
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V3  j\  +  2  (E/N0) 
nBs  y  2(E/N0f 

This  result  holds  for  a  spectrum  that  is  rectangular  over  the  band  Bs. 

Since  the  line  of  constant  DTOA  corresponds  to  an  LOP  on  the  earth,  the  variance  of  the 
DTOA  estimate  can  be  directly  related  to  the  variance  of  the  location  estimate  (along  each  axis) 
by  the  formula  [2] 

[  ^OPTOA  +  ^I  I1'2  (g) 

4  sin2  (0/2)  -  (cos  02  -  cos  0X)2 

where  c  is  the  speed  of  light,  ap  is  the  1-cr  error  in  each  of  the  interceptor  coordinates,  0,  is  the 

pointing  angle  from  the  transmitter  to  receiver  i,  and  0  =  |02  —  0i  |  is  the  angle  separating  the  two 
receivers  when  viewed  from  the  transmitter. 

For  the  case  of  very  accurate  measurements  of  the  receivers’  positions,4  and  the  small  sub¬ 
tending  angle  8  that  is  of  interest  at  high  frequencies.  Equation  (6)  can  be  simplified  by  using 
Up  =  0,  cos 0i  =  cos 02,  and  sin 0/2  ~  0/2: 


COPTOA 

0 


(7) 


Assume  that  the  separation  angle  0  is  equal  to  the  half-power-beamwidth  (HPBW)  of  the 
transmitting  antenna.5  The  HPBW  is  given  by  the  usual  formula 

HPBW  =  1.22  radians  ,  (8) 

where  /  is  the  signal’s  frequency  and  D  is  the  diameter  of  the  circular  aperture  antenna.  Substitut¬ 
ing  this  equation,  along  with  (5),  into  (7),  and  assuming  two  similar,  independent  and  orthogonal6 


4 This  would  be  essentially  true  if  skin-tracking  radar  was  used. 

5The  variation  of  geolocation  performance  with  separation  angle  is  considered  in  Section  5. 

6If  the  LOPs  (labeled  A  and  B ,  say)  are  not  orthogonal,  but  rather  have  an  angle  of  intersection 
0,  then  the  resulting  l-<7  error  is  a  =  /  -J sin  0 . 
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LOPs  due  to  two  DTOA  estimates,  the  final  1  -o  geolocation  error  is 


=  y/2ox  =  2.3  x  10J 


Bs  ■  HPBW  y  2(E/N0) 


1  +  2  (E/No) 


km 


=  6.4x10 ^ 


B. 


2  (E/NoY 


(9) 


where  D  is  given  in  meters,  and  /  and  Ba  are  given  in  hertz.  For  the  sake  of  generality,  the  form 
given  in  (9)  will  be  used  henceforth. 


Similar  equations  can  be  derived  for  the  geolocation  error  based  on  DD  estimates.  In  either 
case,  whether  DTOA-  or  DD-based  data  is  used,  averaging  over  N  independent  sets  of  measurements 
decreases  the  location  error  by  a  factor  of  1  /\fN. 
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3.  STRAWMAN  INTERCEPTOR  SYSTEM 


To  geolocate  a  transmitting  terminal,  at  least  two  independent  LOPs  have  to  be  obtained. 
This  can  be  done  simultaneously  with  only  two  interceptors  (provided  both  DTOA  and  DD  esti¬ 
mates  are  used)  or  with  three  interceptors  (if  only  one  type  of  estimate  is  used).  Using  successive 
gatherings  of  the  DTOA  estimate,  it  is  possible  to  unambiguously  locate  the  terminal  with  only 
two  interceptors,  provided  they  are  moving  with  respect  to  the  terminal.  Indeed,  the  latter  method 
can  be  applied  with  only  one  interceptor  if  the  ephemerides  of  both  the  intended  satellite  and  the 
interceptor  are  accurately  known,  as  are  the  epoch  timing  boundaries  of  the  system.  This  report, 
however,  concentrates  on  the  situation  of  multiple  interceptors  providing  DTOA  estimates.  No 
assumptions  of  accurate  satellite  ephemerides  or  knowledge  of  epoch  timing  are  made,  though  the 
ephemerides  of  the  interceptors  are  assumed  to  be  very  accurately  known. 

Figure  2  illustrates  the  two  specific  interceptor  configurations.  The  first  configuration  consists 
of  four  interceptors  placed  at  the  vertices  of  a  box  (approximately),  centered  on  the  intended 
satellite.  Similarly,  in  the  second  configuration,  three  interceptors  form  an  equilateral  triangle 
about  the  satellite.  In  both  the  box  and  the  triangle  configurations,  which  will  be  symbolized  by  □ 
and  A,  respectively,  the  pair-wise  angular  spacing  of  the  interceptors  is  denoted  by  0.  The  distance 
between  each  interceptor  and  the  satellite  is  0/2  for  the  □  scenario  and  0/\/3  for  the  A  scenario. 

For  each  of  these  configurations  there  is  an  optimal  value  of  6  that  allows  the  most  accurate 
estimation.  As  shown  in  Section  5,  a  value  of  0  =  1  HPBW  is  very  nearly  optimal  for  both  the  □  and 
the  A  scenarios.  Furthermore,  calculations  indicate  that  geolocation  based  on  the  A  configuration 
is  13.5  percent  less  accurate  than  one  based  on  the  □  configuration,  corresponding  to  an  effective 
loss  in  geolocation  SNR  of  1.1  dB. 

The  strawman  interceptor  system  will  be  the  □  configuration  with  1  HPBW  spacing  between 
interceptors.  This  choice  maximizes  the  estimation  accuracy,  though  it  clearly  requires  more  in- 
orbit  assets.  Using  the  formulas  and  graphs  of  Section  5,  it  is  a  straightforward  matter  to  extend 
the  results  from  that  of  the  strawman  system  to  a  □  or  A  configuration  of  arbitrary  interceptor 
spacing. 

Assuming  the  cross-correlation  signed  peak  can  be  detected  in  all  cases,  source  location  ac¬ 
curacy  can  be  enheinced  by  correlating  longer  spans  of  data  eind/or  averaging  successive  DTOA 
estimates.  Frequency  hopping  by  the  source  terminal  seriously  impedes  the  geolocation  procedure 
only  if  the  resulting  SNRs  of  the  hops  seen  at  the  interceptors  are  too  low  to  allow  reliable  detection 
of  the  peak  of  the  ambiguity  function.  In  such  a  case,  it  would  be  necessary  to  collect  and  sort 
multiple  hops  from  all  the  active  sources  in  the  monitored  area,  a  problematical  task.  On  the  otner 
hand,  if  the  SNRs  are  high  enough,  a  DTOA  estimate  can  be  derived  from  each  hop’s  data  and 
then  the  numerous,  albeit  crude,  one-hop  estimates  can  be  averaged  to  yield  the  desired  accuracy. 
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Figure  2.  Strawman  interceptor  system. 
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4.  LINK  BUDGET 


In  the  strawman  interceptor  system,  the  geolocation  accuracy  formula  given  in  Equation  (9)  is 
a  function  of  three  parameters:  the  user  signal  bandwidth  Ba,  the  input  SNR,  and  the  transmitting 
terminal’s  uplink  HPBW.  The  latter  two  parameters  are  quantified  below  for  various  terminals. 

The  received  power-to-noise  ratio  at  a  satellite  from  a  terrestrial  terminal  is 

( Pr /-No ) satellite  =  {EIRP) terminal  ~  lOlog(fcT)  -  Lmisc 

— 201ogil  -  101og(47r)  +  lOlogA  dB-Hz  (10) 

where  k  is  Boltzmann’s  constant  (-228.6  dBW/K-Hz);  T  is  the  satellite  system  temperature;  Lmiac 
is  a  term  into  which  all  the  miscellaneous  losses  are  folded;  R  is  the  path  length;  and  A  is  the 
effective  satellite  antenna  aperture  area.  Assuming  T  =  1000  K,  Lmi3C  =  2  dB,  R  —  40,000  km, 
and  a  couple  of  reasonably-sized  satellite  apertures,  (10)  becomes 


(  Pt  /-No  )  satemte 


(. EIRP)terminal  +  20  dB-Hz,  for  A  =  430cm2 
( EIRP)terminal  +  5  dB-Hz,  for  A  =  15cm2 


(11) 


The  apertures  and  system  temperatures,  which  are  assumed  to  be  the  same  for  both  the 
intended  satellite  and  the  interceptor  satellites,  should  be  easily  attainable. 

Table  1  enumerates  the  EIRPs  and  the  uplink  half-power-beamwidths  of  a  couple  of  proto¬ 
typical  terminals. 


TABLE  1 


Prototypical  Terminal  Characteristics 


Terminal 

Terminal 

Uplink 

Type 

EIRP 

HPBW 

(dBW) 

(radians) 

Type  A 

60 

0.014 

Type  B 

45 

0.062 

Assuming  the  hopping  rate  is  1000  hops-per-second,  Table  2  lists  the  received  SNR  per  hop 
as  a  function  of  the  terminal  type  and  the  effective  satellite  aperture.  Note  that  the  full-power 
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SNR  is  sufficient  to  allow  reliable  peak-detection  based  on  a  single  intercepted  hop’s  worth  of  data, 
mostly  by  a  substantial  margin.  Reduction  of  this  margin  would  entail  the  use  of  lower-power 
transmissions,  which  could  necessitate  the  terminal’s  use  of  lower  data  rates. 


TABLE  2 

Uplink  Signal-to- Noise  Ratios 


Effective 

Uplink  E/N0  (dB) 

Satellite 

Aperture 

Type  A  Terminal 

Type  B  Terminal 

430  cm2 

50 

35 

5  cm2 

35 

20 
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5.  OPTIMUM  ANGULAR  SEPARATION  OF  INTERCEPTORS 


Until  now,  we  have  assumed  the  separation  angle  9  to  be  equal  to  the  transmitting  antenna’s 
HPBW.  This  section  investigates  the  effect  on  the  geolocation  error  a  using  a  different  angle  9. 

The  angle  9  enters  the  geolocation  calculations  directly,  as  an  indicator  of  the  baseline  of 
an  interceptor  pair,  and  indirectly,  through  the  reduction  in  the  received  SNR  caused  by  the 
misdirection  of  the  transmitter’s  antenna  gain  pattern.  These  are  opposing  effects  in  that  an 
increased  9  tends  to  both  enhance  the  estimator  accuracy  (by  increasing  the  baseline)  and  degrade 
it  (by  reducing  the  received  SNR  and  making  it  more  difficult  to  detect  the  peak  of  the  ambiguity 
function).  The  conglomerate  effect  is  quantified  below. 

5.1  Optimum  Angle  for  a  J\(x)/x  Transmitter  Antenna  Pattern 
5.1.1  Box  (□)  Interceptor  Configuration 

First,  consider  the  ideal  power  gain  pattern  of  a  uniformly-illuminated  circular  aperture: 


G(<P)  = 


•J!  (3.234)1 2 

(3.234)  . 


(12) 


where  4  is  the  angular  offset  in  HPBW,  and  the  peak  gain  has  been  normalized  to  unity.  As 
demonstrated  by  the  gains  of  the  first  five  sidelobes  (-17.6,  -23.8,  -28.0,  -31.1,  and  -33.7  dB, 
respectively),  the  gain  falls  off  fairly  rapidly.  Intuitively,  it  may  seem  that  this  quick  decline, 
coupled  with  the  narrow  beamwidths  seen  at  high  frequencies,  forces  the  interceptors  to  be  near 
the  peak  of  the  main  lobe  if  they  are  to  be  at  all  effective.  This  is  indeed  the  case  if  the  signal  is 
undetectable  when  seen  through  a  sidelobe,  thereby  forcing  the  main  lobe  to  be  used.  However, 
if  the  signal  level  is  sufficiently  high  to  allow  reasonable  detection  through  a  sidelobe,  geolocation 
using  that  sidelobe’s  signal  is  possible,  though  it  is  certainly  not  as  efficient  as  using  the  main  lobe’s 
signal. 

As  is  evident  from  Equation  (9),  the  geolocation  variance  a2  is  proportional  to  [92SNR]~  , 
where  SNR  is  an  implicit  function  of  9  by  way  of  the  transmitting  terminal’s  gain  pattern.  Denoting 
the  (power)  gain  pattern  by  G(-),  the  following  proportionality  relation  can  be  written: 


92SNR  a  92  G(9/ 2) 


(13) 


When  the  interceptor  spacing  9  is  changed  from,  say,  0,^  to  9,  the  relative  change  in  the  SNR 
needed  for  an  equivalent  geolocation  variance  is 
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ASNR  = 


6>2  G{6/ 2) 

G(^d/2) 


Using  Odd  =  1  HPBW  as  the  reference  point,  Equation  (14)  becomes 


ASNR  =  2  02  G(0/2)  , 


where  0  is  given  in  units  of  HPBW. 

Figure  3  plots  7  Equation  (15),  as  well  as  G(6/ 2),  as  a  function  of  6. 
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Figure  3.  Relative  geolocation  SNR  as  a  function  of  angular  separation;  J\(x)/x  antenna 
pattern;  □  configuration. 


7Unless  otherwise  stated,  a  J\{x)/x  antenna  pattern  is  assumed. 


This  figure  shows  the  geolocation  efficiency  of  using  various  values  of  9  relative  to  that  of 
using  9  —  1  HPBW,  provided  that  detection  is  possible  at  an  angular  offset  of  9/2  from  the  peak  of 
the  transmitter’s  antenna  pattern  (which  is  pointing  at  the  intended  satellite).  The  latter  condition 
is  necessary  since  an  increased  interceptor  baseline  is  useless  if  the  signal  is  undetectable. 

The  peak  geolocation  efficiency  occurs  at  9  ss  1  HPBW,  while  the  peaks  of  subsequent  ef¬ 
ficiency  lobes  occur  at  about  the  locations  that  correspond  to  the  sidelobes  of  the  gain  pattern. 
That  is,  the  first  four  efficiency  peaks  (after  that  of  the  main  lobe)  are  at  9  =  3.2,  5.2,  7.2,  and 
9.2  HPBW  respectively,  with  corresponding  relative  SNR  degradations  of -4.5,  -6.5,  -7.8,  and  -8  8 
dB.  If  a  sidelobe’s  signal  is  detectable,  then  that  signal  can  also  be  used  for  geolocation,  with  much 
of  the  apparent  loss  in  SNR  (relative  to  the  optimum  of  9  =  1  HPBW)  being  compensated  by  the 
increase  in  the  baseline. 

5.1.2  Triangle  (A)  Interceptor  Configuration 

The  previous  calculations  have  assumed  a  □  interceptor  configuration,  as  defined  in  Section 
3.  Similar  results  can  be  derived  for  the  A  configuration,  using  a  variation  of  Equation  (15)  that 
takes  into  account  the  change  of  the  interceptor-to-satellite  spacing  from  9/2  to  6 /y/3: 

ASNR  =  2.57  92  G(9/VZ)  (16) 

Evaluation  of  this  equation  indicates  the  optimal  efficiency  and  the  associated  9  («  1  HPBW)  to 
be  about  the  same  as  for  the  □  configuration.  For  the  A  array,  this  implies  the  interceptor  SNR 
is  G^l/v^)  =  -4.1  dB  from  the  peak  seen  at  the  intended  satellite,  rather  than  -3  dB  from  the 
peak  for  the  □  array.  Of  course,  the  spacings  between  and  the  values  of  the  subsequent  peaks  are 
different  from  the  □  case. 

Comparing  the  optimal  geolocation  performance  of  the  A  interceptor  array  with  that  of  the 
□  array,  with  9opt  =  1  HPBW  in  either  case, 

*7  A  ,opt 
&a,opt 

This  13.5  percent  degradation  in  geolocation  accuracy  corresponds  to  an  equivalent  loss  in  SNR 
of  1.1  dB,  so  that  ~  29  percent  more  samples  need  to  be  averaged  in  the  A  array  to  equal  the 
performance  of  the  □  array. 
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5.2  Optimum  Angle  for  a  J3(x)/x3  Transmitter  Antenna  Pattern 

It  was  demonstrated  in  Section  5.1  that,  for  a  transmitter  having  a  typical  antenna  gain 
pattern  [such  as  the  ideal  J\{x)/x),  geolocation  using  the  transmitter’s  sidelobes  is  feasible,  provided 
the  SNR  at  each  interceptor  is  high  enough  to  allow  detection  of  the  signal. 

When  the  sidelobes  of  the  transmitting  antenna  gain  pattern  are  reduced,  e.g.,  by  modifying 
the  feed,  the  resulting  decrease  in  detectability  may  be  sufficient  to  restrict  all  practical  interceptor 
separation  angles  to  be  in  the  region  of  the  main  lobe.  In  such  a  case,  realistic  limitations  on  the 
minimum  interceptor- to-satellite  orbital  distance  would  allow  geolocation  of  only  those  terminals 
that  have  small  antennas.  The  terminals  having  larger  antennas  (and  smaller  HPBWs)  would 
require  unrealistically  small  interceptor- to-S/C  orbital  distances. 

As  an  example,  parabolic-squared  tapering  [3]  of  the  antenna  feed  produces  a  low-sidelobe 
gain  pattern  that  has  the  form  J3(x)/x3.  The  first  sidelobe  of  the  J^xj/x'6  pattern  is  more  than 
30  dB  down  from  the  peak,  which  is  about  the  same  level  as  for  the  fourth  sidelobe  of  the  J\  (x) /x 
pattern.  Sidelobe  geolocation  that  is  feasible  by  a  margin  of  less  than  12  dB  for  a  J\(x)/x  pattern 
is  not  feasible  for  the  case  of  a  Jz(x)/x*  pattern.  Mainlobe  geolocation  is  similar  for  both  patterns. 

As  in  the  J\(x)/x  case,  the  relative  degradation  in  the  J^(x)/x3  geolocation  accuracy,  due  to 
a  change  from  a  □  configuration  to  a  A  configuration,  is  about  1  dB,  requiring  about  26  percent 
more  data  for  equivalent  performance. 

5.3  Optimum  Interceptor  Separation 

The  generic  interceptor-separation  results  of  Section  5.1  can  be  applied  to  the  terminals  listed 
in  Section  4.  With  uniformly-illuminated  transmitter  apertures,  a  □  interceptor  configuration,  and 
a  transmitter-to-satellite  distance  of6  R  =  40,000  km,  the  relative  antenna  gains  of  the  two  terminal 
types  are  shown  in  Figure  4  as  a  function  of  the  interceptor-to-S/C  distance  (=s  RO/2).  Figure  5 
graphs  the  relative  geolocation  SNRs  for  the  two  apertures  also  as  a  function  of  the  interceptor-to- 
S/C  distance,  hereafter  referred  to  as  V.  That  is,  in  comparison  with  the  point  T>g- 1  hpbw >  the 
geolocation  accuracy  is  reduced  by  an  amount  equivalent  to  the  drop  in  SNR  depicted  in  Figure  5. 


5.3.1  Optimization  for  Each  Terminal 

Figure  5  illustrates  that  the  optimum  values  of  V  (corresponding  to  a  6  of  1  HPBW)  are 
270  km  and  1230  km  for  the  Type- A  and  Type-B  apertures,  respectively.  These  are  the  individ¬ 
ually  optimum  values  of  P,  assuming  no  restriction  on  the  minimum  distance.  The  imposition 


8A  20°  elevation  look  angle  and  geostationary-altitude  satellites  are  assumed. 
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Figure  4-  Relative  gain  from  various  transmitter  apertures. 


of  a  restriction  on  V  modifies  the  corresponding  values  of  I>opt,  as  demonstrated  by  the  following 
examples. 

If  V  is  restricted  to  be  >  ~  500  km,  V  =  Vhpbw  for  the  Type-B  antenna  as  before,  the 
Type-A  terminal  must  be  detected  through  its  first  sidelobe  (V  =  860  km).  The  effective  loss  in 
geolocation  SNR  of  the  latter  terminal  is  only  4.5  dB,  provided  the  ambiguity  peak  can  be  detected 
despite  the  14.6  dB  reduction  in  gain  relative  to  the  3-dB  pom.' . 

In  a  similar  manner,  extending  the  lower  limit  of  V  to  ~  1000  km  forces  the  Type-A  terminal 
to  be  received  through  its  second  sidelobe.  It  is  essential  to  keep  in  mind  that  the  small  relative 
degradations  in  geolocation  SNR  illustrated  in  Figure  5  assume  that  the  peak  of  the  ambiguity 
function  is  reliably  detectable.  The  practicality  of  detecting  the  peak  is  determined  by  the  curves 
of  Figure  4  coupled  with  the  EIRPs  and  link  budgets  of  Section  4 

5.3.2  Joint  Optimization 

In  the  previous  section,  the  optimum  interceptor-S/C  distance  V  was  considered  separately 
for  each  source  with  a  given  antenna  size,  and  under  various  limitations  on  T>.  To  detect  both 
source  types  simultaneously,  the  jointly  optimum  distance  is  equal  to  the  optimum  distance  of  the 
smaller  aperture,  provided  there  is  sufficient  margin  to  detect  the  ambiguity  peak  on  the  appropriate 
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INTERCEPTOR  ■  S/C  SEPARATION  (km) 

Figure  5.  Relative  geolocation  SNR  for  various  transmitter  apertures. 


antenna  sidelobes.  However,  insufficient  sidelobe  margins  may  force  the  use  of  a  smaller  V  to  allow 
peak  detection,  with  a  concomitant  increase  in  the  amount  of  data  required  to  achieve  the  desired 
geolocation  accuracy. 

Consider  the  joint  geolocation  of  Type-A  and  Type-B  terminals.  Selecting  a  V  of  860  km 
puts  the  interceptor  at  the  peak  of  the  larger  antenna’s  first  sidelobe  and  well  within  the  main  lobe 
of  the  smaller  antenna.  The  relative  antenna  gains  are  -14.6  dB  and  +1.6  dB,  respectively,  and 
the  corresponding  relative  geolocation  SNRs  are  -4.5  dB  and  -1.5  dB.9 

On  the  other  hand,  having  V  —  270  km  lets  the  interceptor  be  at  Vhpbw  for  the  Type-A 
source;  this  distance  places  the  interceptor  near  the  peak  of  the  Type-B  source.  The  relative  gain 
and  geolocation  SNR  of  the  Type-A  source  are  both  0  dB,  while  those  for  the  Type-B  source  are 
+2.9  dB  and  -10.3  dB,  respectively.  In  other  words,  for  the  case  of  the  smaller-aperture  source  the 


9The  gains  referred  to  tire  the  antenna  power  gains  at  a  given  V  relative  to  the  gains  at  Vhpbw, 
or  3  dB  down  from  the  mainlobe  peak.  Similarly,  the  geolocation  SNRs  are  presented  relative  to 
the  optimum  values,  which  occur  at  V  =  Vhpbw- 
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relative  gain  is  increased,  but  the  reduced  baseline  causes  a  reduction  in  the  relative  geolocation 
accuracy,  thereby  extending  the  required  data-gathering. 

The  choice  between  the  above  two  options  for  Q  depends  upon  the  uplink  margins  of  the 
sources,  their  relative  importance,  and  the  limitations,  if  any,  on  the  maximum  assured  duration  of 
data.  The  use  of  these  guidelines  is  illustrated  by  the  following  special  case:  if  the  smaller  sources 
are  more  important,  there  are  no  severe  time  limitations,  and  the  larger  sources  have  uplink  margins 
(to  the  intended  satellite)  of  at  least  ~  20  dB,  then  the  appropriate  choice  for  V  is10  860  km. 


10This  value  is  for  geostationary  satellites  and  interceptors  with  the  appropriate  angular  spacing, 
viewed  from  an  elevation  angle  of  20°.  The  distance  decreases  with  increasing  elevation  angle, 
reaching  780  km  at  90°. 
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6.  GEOLOCATION  OF  TERMINALS:  IDEAL  PERFORMANCE 


The  Cramer-Rao  bounds  on  the  geolocation  error,  detailed  in  Section  2,  will  now  be  used 
to  calculate  the  geolocation  regions  for  the  terminals  detailed  in  Tables  2  and  1.  Throughout 
the  following  it  is  assumed  that  the  □  interceptor  configuration  is  used,  with  an  angular  spacing 
between  opposing  interceptors  of  6  ~  1  HPBW.  The  interceptor- spacecraft  distances  (2?)  associated 
with  this  optimum  angular  spacing  for  Type-A  and  Type-B  aperture  sources  are  at  270  and  1230 
km,  respectively.  The  effect  of  other  nonoptimal  values  of  D  on  the  results  presented  below’  can  be 
calculated  using  the  graphs  previously  given  in  Section  5. 

For  each  terminal  type  and  each  filter  bandwidth  ( Ba ),  Equation  (9)  was  used  to  calculate 
the  1— <r  geolocation  error  for  DTOA-only  measurements  as  a  function  of  the  input  SNR  at  the 
interceptors’  front  ends.  The  geolocation  results,  presented  in  Figure  6,  are  for  the  reception  of  a 
single  hop  having  an  SNR  of  E/No.  As  always,  a  prerequisite  for  this  geolocation  accuracy  is  the 
detection  of  the  peak  of  the  ambiguity  function. 
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Figure  6.  l-o  geolocation  error. 
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When  the  SNR  at  the  interceptor  is  greater  than  about  15  dB,  a  single  hop  is  sufficient  to 
detect  the  peak.  Often,  a  10  dB  SNR  is  used  as  a  rule  of  thumb  to  indicate  reliable  detection  based 
on  a  single  measurement  [lj.  In  a  very  wide  hopping  band,  however,  the  number  of  time-frequency 
cells  that  are  to  be  searched  for  a  detection  is  very  large,  and  the  required  false-alarm  probability 
per  cell  is  extremely  small,  thereby  requiring  a  higher  SNR  (such  as  15  dB)  in  order  to  make  a 
detection  decision  on  the  basis  of  only  a  single  hop.  Nevertheless,  peak  detection  is  possible  at  a 
much  lower  SNR,  but  the  decision  is  then  made  using  multiple  hop  measurements,  and  the  peak 
estimation  is  done  in  a  decision-directed  manner  using  the  same  stored  measurements. 


20 


7.  OBSERVATIONS 


7.1  Interceptor  G/T 

It  was  assumed  in  Section  4  that  the  receiver  antenna  apertures  of  the  intended  satellite 
and  each  interceptor  were  rather  small.  However,  the  type  of  antenna  employed  by  a  practical 
interceptor  could  range  from  a  small-aperture  earth-coverage  horn  to  a  larger-aperture  narrow 
spot  beam  antenna.  If  the  interceptor’s  antenna  aperture  was  considerably  greater  than  that  of 
the  intended  satellite,  the  intercept  region  of  the  terminals  transmitting  into  that  beam  would  be 
substantially  increased. 

For  example,  if  a  terminal  is  transmitting  both  into  the  intended  satellite’s  earth-coverage 
beam  and  into  an  interceptor’s  spot  beam,  the  transmitter  would  be  very  easily  detected  and 
geolocated.  The  interceptor’s  considerable  detection  margin  could  be  used  either  to  geolocate 
terminals  through  the  sidelobes  of  their  uplink  beams  or  to  restrict  the  search  window  in  the 
manner  indicated  in  the  following  section. 

An  additional  benefit  to  the  interceptors  using  highly  directional  antennas  is  the  inherent 
spatial  discrimination  that  they  would  provide,  making  it  easier  to  separate  the  presence  of  a 
number  of  transmitters  in  the  region  being  monitored. 

On  the  other  hand,  if  the  interceptor  satellites  are  using  broad  uplink  beams  while  the  intended 
satellite  has  a  pencil  beam,  the  signals  received  by  the  interceptors  will  be  much  weaker,  and 
geolocation  will  be  more  difficult. 

7.2  Restricting  the  Search  Window 

In  many  of  the  likely  scenarios  the  SNR  received  by  the  interceptors  will  be  substantial,  i.e., 
greater  than  15  dB,  and  the  number  of  hops  needed  for  geolocation  will  not  be  very  large.  Therefore, 
the  search  procedure  can  be  considerably  restricted  in  any  of  several  dimensions  (namely,  frequency, 
space,  and  time)  and  still  provide  very  quick  geolocation. 

First,  and  probably  most  importantly,  only  a  portion  of  the  hopping  band  may  be  monitored. 
If  the  transmitting  terminal  hops  uniformly  across  the  frequency  band,  enough  data  can  si.*l  be 
collected  to  form  a  good  geolocation  estimate. 

Second,  in  the  case  of  the  interceptors  having  highly  directive  beams,  they  can  be  coordinated 
to  scan  across  a  region.  Because  for  geolocation  at  moderate-to-high  SNR  only  a  few  hops  are 
necessary,  the  periodic  coverage  of  the  entire  area  would  be  adequate  to  catch  transmissions  that 
last  more  than  a  few  seconds. 

Finally,  the  processing  in  the  ground  station  need  not  be  done  in  real  time.  The  station  may 
take,  say,  1  s  to  process  100  ms  worth  of  data,  for  a  duty  cycle  of  9  percent.  In  many  cases,  this 
partial  temporal  coverage  will  provide  adequate  geolocation  service. 


8.  CONCLUSION 


The  results  of  the  theoretical  calculations  detailed  in  this  report  indicate  that  geolocation 
of  frequency-hopping  terminals  by  a  network  of  spaceborne  interceptors  and  an  accompanying 
terrestrial  processing  station  is  feasible  under  the  indicated  conditions.  The  optimum  angular 
separation  of  the  interceptor  satellites  was  found  to  be  approximately  equal  to  the  uplink  HPBW 
of  the  transmitting  antenna,  and  the  effect  of  nonoptimal  separation  was  discussed.  A  system  to 
geolocate  several  terminal  types  would  require  a  compromise  choice  for  the  interceptor  separation 
because  different  types  of  terminals  have  different  uplink  HPBWs.  The  extent  of  the  compromise 
would  depend  on  the  relative  value  afforded  by  the  interception  of  the  various  terminals. 

As  a  final  point,  it  should  be  noted  that  considerable  effort  is  involved  in  forming  the  geolo¬ 
cation  system  type  detailed  in  this  report.  To  create  such  a  system  it  would  be  necessary  to: 

•  Place,  maintain  in  appropriate  orbit,  and  coordinate  several  intercepting  transpon¬ 
ders  near  each  intended  satellite,  thereby  requiring  the  expenditure  of  substantial 
resources  to  form  complete  interception  coverage  of  a  numerous-satellite  system. 

•  Ensure  that  the  interceptors’  G/T  are  not  much  less  than  the  G/T  of  the  intended 
satellite  uplink. 

Even  for  an  ideal  system  of  this  kind,  there  may  be  situations  where  the  desired  geolocation  accuracy 
and  processing  speed  are  not  readily  available.  For  example,  it  would  be  difficult  to  distinguish 
multiple  terminals  that  are  simultaneously  transmitting  into  a  high-gain  uplink  beam  using  a  low- 
data-rate  mode  near  its  threshold.  In  this  case,  numerous  hops  would  be  needed  to  find  the  peak  of 
each  terminal’s  ambiguity  function,  and  much  raw  data  would  have  to  be  stored  before  the  various 
terminals  could  be  differentiated. 
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